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ABSTRACT 

The Wilkinson Microwave Anisotropy Probe (WMAP) has mapped the full sky in Stokes /, Q, and U parameters 
at frequencies 23, 33, 41, 61, and 94 GHz. We detect correlations between the temperature and polarization maps 
significant at more than 10 standard deviations. The correlations are inconsistent with instrument noise and are 
significantly larger than the upper limits established for potential systematic errors. The correlations are present in 
all WMAP frequency bands with similar amplitude from 23 to 94 GHz, and are consistent with a superposition of 
a CMB signal with a weak foreground. The fitted CMB component is robust against different data combinations 
and fitting techniques. On small angular scales (0 < 5°), the WMAP data show the temperature-polarization 
correlation expected from adiabatic perturbations in the temperature power spectrum. The data for t > 20 agree 
well with the signal predicted solely from the temperature power spectra, with no additional free parameters. We 
detect excess power on large angular scales (0 > 10°) compared to predictions based on the temperature power 
spectra alone. The excess power is well described by reionization at redshift 11 < z r < 30 at 95% confidence, 
depending on the ionization history. A model-independent fit to reionization optical depth yields results consistent 
with the best-fit ACDM model, with best fit value r = 0.17 ± 0.04 at 68% confidence, including systematic and 
foreground uncertainties. This value is larger than expected given the detection of a Gunn-Peterson trough in the 
absorption spectra of distant quasars, and implies that the universe has a complex ionization history: WMAPhas 
detected the signal from an early epoch of reionization. 


1. INTRODUCTION 

Line.ar polarization of the cosmic microwave background 
(CMB) results from anisotropic Thomson scattering of CMB 
photons by free electrons. By symmetry, an isotropic radiation 
field can not generate a net polarization. Any net polarization 
results from the quadrupole moment of the CMB temperature 
distribution seen by each scattered Multiple scattering sup- 
presses polarization by damping the temperature anisotropy; 
hence, CMB polarization originates primarily from epochs 
when the opacity was of order unity or less. Standard cos- 
mological models predict two such epochs, corresponding to 
two characteristic angular scales. The first is the decoupling 
surface at redshift z ~ 1088, when the ionization fraction x e 
abruptly falls from near unity to near zero. The acoustic horizon 
at decoupling subtends an angle 6 « 1°; polarization on these 
scales reflects conditions in the photon-baryon fluid just prior 
to recombination. Polarization data from decoupling comple- 
ment measurements of the temperature anisotropy. Astrophys- 
ical sources generate additional polarization as ionizing radia- 
tion from the first collapsed objects reionizes the intergalactic 


medium. For reionization at redshift z < 50 the horizon is on 
large angular scales, 0 > 5°. Polarization on these scales di- 
rectly probes the poorly-understood process of reionization. 

Since CMB polarization originates at modest opacity, the 
underlying temperature anisotropy is not heavily damped and 
remains observable today. Precise predictions can be made 
of the average polarization pattern expected from a given 
power spectrum of temperature anisotropy (Rees (1968); Kaiser 
(1983); Bond & Efstathiou (1984); Coulson et al. (1994); 
Kamionkowski et ai. (1997); Zaldarriaga & Seljak (1997); Hu 
& White (1997); for recent reviews, see Kosowsky (1996); Hu 
& Dodelson (2002)). The pattern of polarization on the sky is a 
vector field with both an amplitude and direction at each point, 
and can be separated into two scalar fields, one giving the curl 
and the other the gradient component (called B and E modes in 
analogy with electromagnetic fields). The DASI collaboration 
has detected CMB polarization on angular scales ~ 0.°5 (Ko- 
vac et al. 2002). DASI reports an E mode signal significant at 
4.9a and a TE temperature-polarization correlation significant 
at 2a. Both signals are consistent with the “concordance” cos- 
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mological model (spatially flat model dominated by a cosmo- 
logical constant and cold dark matter; see, e.g„ Hu & Dodelson 
(2002)) and support an adiabatic origin for the CMB tempera- 
ture anisotropy. 

The Microwave Anisotropy Probe has mapped the full sky in 
the Stokes /, Q , and U parameters on angular scales 6 > 0.°2 
in 5 frequency bands centered at 23, 33, 41, 61, and 94 GHz 
(Bennett et al. 2003a). WMAP was not designed solely as a 
polarimeter, in the sense that none of its detectors are sensitive 
only to polarization. Incident radiation in each differencing as- 
sembly (DA) is split by an orthomode transducer (OMT) into 
two orthogonal linear polarizations (Page et al. 2003b; Jarosik 
et al. 2003). Each OMT is oriented so that the electric field 
directions accepted in the output rectangular waveguides lie at 
±45° with respect to th eyz symmetry plane of the satellite (see 
Bennett et al. (2003a) Fig. 2 for the definition of the satellite 
coordinate system). The two orthogonal polarizations from the 
OMT are measured by two independent radiometers. Each ra- 
diometer differences the signal in the accepted polarization be- 
tween two positions on the sky (the A and B beams), separated 
by ~ 140°. 

The signal from the sky in each direction h can be decom- 
posed into the Stokes parameters 

T C h ) = 7(n) + Q(h) cos I'y + U (n) sin 27, ( 1 ) 

where we define the angle 7 from a meridian through the Galac- 
tic poles to the projection on the sky of the E-plane of each 
output port of the OMT (Fig. 1). In principle, by tracking the 
orientation of the OMTs on the sky as the satellite scan pat- 
tern observes each sky pixel in different orientations, each ra- 
diometer could independently produce a map of the Stokes /, 
Q, and U parameters. In practice, the non-uniform coverage of 
7 at each pixel would generate significant correlations between 
the fitted Stokes parameters, allowing leakage of the dominant 
temperature anisotropy into the much fainter polarization maps. 
We avoid this problem by differencing the outputs of the two 
radiometers in each differencing assembly in the time-ordered 
data. Denoting the two radiometers by subscripts 1 and 2, the 
instantaneous outputs are 

ATi = I(h A ) + Q(,n A )cos2'y A + U(n A )sin2'y A 

- ~Q(Ab) COS 27 B - U (n B ) sin 2^ B (2) 

and 

AT 2 = I(h A )- Q(h A ) cos 2~{ A - U (n A ) sin 274 

- I(fi B )+ Q(n B )cos2'y B +U(fi B )sin2y B . 

The sum 

AT, = l -(AT + AT 2 ) = I(n A )-I(h B ) (3) 

is thus proportional to the unpolarized intensity, while the dif- 
ference 

AT P = l -(AT l -AT 1 ) 

= QriP-A ) cos 274 + u (h A ) sin 27^ 

- Q(kb) cos 2jb — U (kb) sin 275. (4) 

is proportional only to the polarization. We produce full-sky 
maps of the Stokes /, Q , and U parameters from the sum and 
difference time-ordered data using an iterative mapping algo- 
rithm. Since the polarization is faint, the Q and U maps are 
dominated by instrument noise and converge rapidly (Hinshaw 
et al. 2003a). 


The Stokes Q and U components depend on a specific choice 
of coordinate system. For each pair of pixels, we define 
coordinate-independent quantities 

Q'=Qcos(20)-ff/sin(20) 

U' = U cos(2 (f>) - Q sin(2<£), (5) 

where the angle 0 rotates the coordinate system about the 
outward-directed normal vector to put the meridian along 
the great circle connecting the two positions on the sky 
(Kamionkowski et al. 1997; Zaldarriaga & Seljak 1997). All 
of our analyses use these coordinate-independent linear combi- 
nations of the Q and U sky maps. 

Simulations of the mapping algorithm demonstrate that 
WMAP can accurately recover the polarization pattern on the 
sky, even after allowing for residual calibration uncertainty in 
the individual radiometer channels. However, non-ideal instru- 
mental signals affect the Q and U sky maps to a greater extent 
than the unpolarized / maps. The spacecraft spin about its z 
axis sweeps the beams across the sky in a direction 45° from 
the OMT orientation, preferentially coupling signals not fixed 
on the sky into the U map. Residual striping exists to a lesser 
extent in the I and Q maps. Systematic errors in the individual 
Q and U maps are not yet fully assessed; consequently, we de- 
fer detailed analysis of the Q or U maps to a later paper. Cross- 
correlations between maps are largely unaffected by striping or 
any other channel-specific signal, allowing much simpler analy- 
sis of the faint polarization signal than would be possible for the 
individual QotU maps. This paper discusses the temperature- 
polarization (TE) correlation in the WMAP one-year sky maps. 

We compute the temperature-polarization cross-correlation 
using three different techniques: the two-point correlation func- 
tion, a quadratic estimator for the power spectrum, and a “tem- 
plate” comparison in pixel space between the polarization maps 
and the predicted polarization given the observed pattern of 
temperature anisotropy. All three methods yield similar results 
despite disparate treatments of the data. 


2. CORRELATION FUNCTION 


The simplest measure of temperature-polarization cross- 
correlation is the two-point angular correlation function 


C ,q {9) = 


T,ij I i&j w i w j 


( 6 ) 


where / and j are pixel indices and w are the weights. To avoid 
residual effects of 1 // noise, we force the temperature map 
to come from a different frequency band than the polarization 
maps, and thus use the temperature map at 61 GHz (V band) 
for all correlations except the V-band polarization maps, which 
we correlate against the 41 GHz (Q band) temperature map. 
Since WMAP has a high signal-to-noise ratio measurement of 
the CMB temperature anisotropy, we use unit weight (vv/ = 1) 
for the temperature maps and noise weight (wj = Nj/afy for the 
polarization maps, where Nj is the effective number of obser- 
vations in each pixel j and < 7 0 is the standard deviation of the 
white noise in the time-ordered data (Table 1 of Bennett et al. 
(2003b)). We compare the correlation functions to Monte Carlo 
simulations of a null model, which simulates the temperature 
anisotropy using the best-fit ACDM model (Spergel et al. 2003) 
but forces the polarization signal to zero. Each realization gen- 
erates a CMB sky in Stokes /, Q, and U parameters, convolves 
this simulated sky with the beam pattern for each differencing 
assembly, then adds uncorrelated instrument noise to each pixel 
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Fig. I . — Geometry for Stokes Q and U parameters. WMAP measures polarization by differencing two orthogonal polarization channels, then solving for Q and 
U as the spacecraft compound spin projects the OMT onto the sky at different angles 7 relative to the Galactic meridians. All analysis uses coordinate-independent 
quantities 0 and U' defi ned with respect to the great circle connecting a pair of pixels (see text). 


in each map. We then co-add the simulated skies in each fre- 
quency band and compute C IQ (0) using the same software for 
both the WMAP data and the simulations. All analysis uses 
only pixels outside the WMAP KpO foreground emission mask 
(Bennett et al. 2003c), approximately 76% of the full sky. 

Figure 2 shows C IQ (0) derived by co-adding the individual 
correlation functions for the frequencies 41, 61, and 94 GHz 
(Q, V, and W bands) least likely to be affected by Galactic 
foregrounds. The grey band shows the 68 % confidence inter- 
val for the null simulations. It is clear that WMAP detects a 
temperature-polarization signal at high statistical confidence, 
and that signals exist on both large and small angular scales. 
We define a goodness-of-fit statistic 

x 2 =£ tC P -(c s ? m >LM-' A-(d s i)] b , (?) 

ab 

where C^ AP is the co-added correlation function from WMAP 
data, (C^) is the mean from the Monte Carlo simulations, 
and M is the covariance matrix between angular bins a and 
b derived from the simulations. We find x 2 - 207 for 78 de- 
grees of freedom when comparing WMAP to the null model: 
WMAP detects temperature-polarization correlations signifi- 
cant at more than 10 standard deviations. 

2.1. Systematic Error Analysis 

Having detected a significant signal in the data, we must de- 
termine whether this signal has a cosmological origin or re- 
sults from systematic errors or foreground sources. We test the 
converg ence of the mapping algorithm using end-to-end simu- 
lations, comparing maps derived from simulated time-ordered 
data to the input maps used to generate the simulated time se- 
ries. The simulations include all major instrumental effects, 
including beam ellipticity, radiometer performance, and instru- 
ment noise (including 1 //component), and are processed using 
the same map-making software as the WMAP data (Hinshaw 


et al. 2003a). The Q and U maps converge rapidly, within the 
30 iterations required to derive the calibration solution. Cor- 
relations in the time-ordered data introduce an anti-correlation 
in the U map at angles corresponding to the beam separation, 
with amplitude 0.5% of the noise in the map. This effect is in- 
dependent for each radiometer and does not affect temperature- 
polarization cross-correlations. Similarly, residual 1 // noise in 
the time series can create faint striping in the maps, but does 
not affect cross-correlations. 

The largest potential systematic error in the temperature- 
polarization cross-correlation results from bandpass mis- 
matches in the amplification/detection chains. We calibrate the 
WMAP data in thermodynamic temperature using the Doppler 
dipole from the satellite’s orbit about the Sun as a beam- 
filling calibration source (Hinshaw et al. 2003a). Astrophysi- 
cal sources with a spectrum other than a 2.7 K blackbody are 
thus slightly mis-calibrated. The amplitude is dependent on the 
product of the source spectrum with the unique bandpass of 
each radiometer. If the bandpasses in each radiometer were 
identical, the effect would cancel for any frequency spectrum, 
but differences in the bandpasses between the two radiometers 
in each DA generate a non-zero residual in the difference signal 
used to generate polarization maps (Eq. 4). This signal is spa- 
tially correlated with the unpolarized foreground intensity but 
is independent of the orientation of the radiometers on the sky 
(polarization angle 7 ). In the limit of uniform sampling of 7 this 
term drops out of the sky map solution. However, the WMAP 
scan pattern does not view each pixel in all orientations; unpo- 
larized emission with a non-CMB spectrum can thus be aliased 
into polarization if the bandpasses of the two radiometers in 
each DA are not identical. This is a significant problem only at 
23 GHz (K band), where the foregrounds are brightest and the 
bandpass mismatch is largest. 

We quantify the effect of bandpass mismatch using end-to- 
end simulations. For each time-ordered sample, we compute 
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Fig. 2. — Temperature-polarization correlation function for WMAP co-added QVW data. The gray band shows the 68% confi dence interval for similar co-added 
data taken from Monte Carlo simulations without polarization. The inset shows data for $ < 10°. The data are inconsistent with no temperature-polarization 
cross-correlations at more than 10 standard deviations. Note that the data are not independent between angular bins. 


the signal in each radiometer using an unpolarized foreground 
model and the measured pass bands in each output channel 
(Jarosik et al. 2003). We then generate maps from the simu- 
lated data using the WMAP one-year sky coverage and com- 
pute C 1Q {6) using the output /, Q , and U maps from the simu- 
lation. Figure 3 shows the predicted signal at K band. We treat 
this as an angular template and compute the least-squares fit of 
the WMAP data to this bandpass template to determine the am- 
plitude of the effect in the observed correlation functions. We 
correct the WMAP correlation functions C IQ {6) and C IU (0) at 
K and Ka bands by subtracting the best-fit template amplitudes. 
The fitted signal has peak amplitude of 8 fiK 2 at 23 GHz and 5 
fiK 2 at 33 GHz. No other channel has a statistically significant 
detection of this effect. 

Sidelobe pickup of polarized emission from the Galactic 


plane can also produce spurious polarization at high latitudes 
in the Q and U maps. We estimate this effect using the mea- 
sured far-sidelobe response for each beam in each polarization 
(Barnes et al. 2003). The simplest approach would be to es- 
timate the signal in each time-ordered sample, convolving the 
full sky sidelobe response with the Stokes /, Q , and U maps 
given the instantaneous orientation of the beams for each sam- 
ple. Such an approach is computationally expensive. We in- 
stead approximate the signal in each pixel by convolving the 
full sky sidelobe response with the one-year Q and U maps. 
For each pixel, we fix one beam on that pixel while sweeping 
the other beam through all orientations achieved in flight. The 
average from the convolution yields the sidelobe contribution 
for the pixel in question. Details of the sidelobe maps are pre- 
sented in Barnes et al. (2003). We correlate the sidelobe maps 


Table 1 

Null Tests for IQ and IU Sum and Difference Data 2 


Correlation 

Range 

DOF 

— Sum Map b — 
X 2 Probability d 

— Difference Map c — 
X 2 Probability d 

IQ 

9 <5° 

20 

62.1 

3x10"* 

23.6 

0.26 

IQ 

9 >5° 

58 

145.1 

2xl0' 9 

66.0 

0.22 

IU 

9 <5° 

20 

30.9 

0.06 

10.8 

0.95 

IU 

9 >5° 

58 

66.1 

* 0.22 

50.4 

0.95 


a X 2 comparison of the WMAP correlation functions C ! ^(9) and C ru ($) to a null hypothesis of CMB temperature anisotropy and instrument noise, but no polarization. 
Temperature-polarization signals of cosmic origin should contribute to only C IQ {$) in the sum maps. All other tests are consistent with the null hypothesis. 
^Polarization sum maps (Ql+Q2)/2, (Vl+V2)/2, (Wl+W2)/2, and (W3+W4)/2 co-added with noise weights. 

^Polarization difference maps (Ql-Q2)/2, (Vl-V2)/2, (Wl-W2)/2, and (W3-W4)/2 co-added with noise weights. 

^Probability to randomly obtain x 2 larger than the measured value. 
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Fig. 3. — Angular templates for potential systematic errors caused by bandpass mismatch between the 2 radiometers in each differencing assembly. We fi t this 
template to the correlation functions from each DA to detect or limit systematic errors related to bandpass mismatch in the main beam. The effect is signifi cant only 
in K and Ka bands, which have the brightest unpolarized foregrounds. 


with the temperature anisotropy maps in each channel to esti- 
mate the systematic error in the temperature-polarization cor- 
relation. Sidelobe pickup of polarized structure in the Galactic 
plane is less than 1 pY? in C IQ (0) at 23 GHz and below 0.1 
fiK 2 in all other bands. The effect of bandpass mismatch in 
the far sidelobes (as opposed to the main beam) is similarly 
weak, with limits 1.3 p K 2 at 23 GHz and less than 0.05 pK 2 
in all other bands. We do not apply any corrections for these 
effects, but simply propagate the associated systematic uncer- 
tainty throughout our analysis. Note that all of these systematic 
errors depend on the Galactic foregrounds, and have different 
frequency dependence than CMB polarization. 

Null tests provide an additional check for systematic errors. 
Thomson scattering of scalar temperature anisotropy produces 
a curl-free polarization pattern. A non-zero cosmological sig- 
nal is thus expected only for the IQ (TE) correlation, whereas 
systematic errors or foreground sources can affect both the IQ 
and IU (TB) correlations. We also test linear combinations of 
radiometer maps which cancel the polarization signal but which 
test for systematic effects. We compute the IQ and IU cor- 
relation functions by correlating the Stokes / sum map from 
the Q- or V-band (as noted above) with the polarization dif- 
ference maps (gl-22)/2, (Vl-V2)/2, (Wl~W2)/2, and 
(W3 -W4)/2. We then co-add the results with their noise 
weights, and compare the co-added result for the polarization 
difference maps to a similar computation for the polarization 
sum maps. The temperature (Stokes I) map in all cases is a sum 
map; the test is thus primarily sensitive to systematic errors in 
the polarization data. 

Table 1 shows results of the null tests. We compare C IQ (0) 
and C IU ( 6 ) for the sum and difference maps to a null hypothesis 
that the data consist of Stokes / and instrument noise, with no 


polarization in the Stokes QorU maps. We break the data into 
2 angular regimes to differentiate between signals at decoupling 
vs reionization. We find a clear signal detection for C lQ {6) in 
the sum map for both angular scales. All other tests are con- 
sistent with instrument noise - there is no evidence for addi- 
tional systematic errors in the temperature-polarization cross- 
correlation. 

2.2. Foregrounds 

Galactic emission is not a strong contaminant for CMB 
temperature anisotropy, but could be significant in polariza- 
tion. WMAP measurements of unpolarized foreground emis- 
sion show synchrotron, free-free, and thermal dust emission 
all sharing significant spatial structure (Bennett et al. 2003c). 
Of these components, only synchrotron emission is expected to 
generate significant polarization; other polarized sources such 
as spinning dust are limited to less than 5% of the total intensity 
at 33 GHz. 

Synchrotron emission from electrons accelerated in the 
Galactic magnetic field is the dominant unpolarized foreground 
at frequencies below ~50 GHz. Although it is known to be lin- 
early polarized, previous radio surveys provide little guidance 
for the high-latitude polarization at mm wavelengths. Extrap- 
olation of radio polarization maps (Brouw & Spoelstra 1976) 
to millimeter wavelengths indicate a polarization fraction be- 
tween 10% and 50% depending on Galactic latitude (Lubin & 
Smoot 1981). The unpolarized component has a steep angu- 
lar dependence, with power spectrum a oc t 1 (Bennett et al. 
2003c). The angular dependence of the polarized component is 
expected to be even steeper (Baceigalupi et al. 2001; Bruscoli 
et al. 2002; Tucci et al. 2002), suggesting that foreground po- 
larization is most likely to affect results on large angular scales. 
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Fig. 4. — Fitted CMB (left) and foreground (right) components from a multi- frequency decomposition of the measured two-point correlation functions. Top panels 
show the IQ (TE) correlation, while bottom panels show IU (TB). The CMB component is shown in units of thermodynamic temperature, while the foreground is 
shown in antenna temperature evaluated at 41 GHz. The best-fi t spectral index is $ = -2.7 ± 0.3. The fi tted CMB and foreground components are stable as additional 
low-frequency data are included in the fi t. 


Radio maps at low Galactic latitude, however, demonstrate 
that the polarization intensity is not necessarily well correlated 
with the unpolarized intensity, complicating template analysis 
for temperature-polarization cross-correlations (Uyamker et al. 
1998, 1999). We thus use the frequency dependence of the 
measured temperature-polarization cross-correlation to sepa- 
rate cosmic from foreground signals. 

Foreground polarization above 40 GHz is faint: fitting the 
correlation functions at 41, 61, and 94 GHz (Q, V, and W 
bands) to a single power-law C 1Q (9,v) = d®(9) (v/vo) 20 yields 
spectral index fi = -0.4 ± 0.4, consistent with a CMB signal 
fi = 0) and inconsistent with the spectral indices expected for 
synchrotron fi « -3), spinning dust fi % -2), or thermal dust 
fi « 2). The measured signal can not be produced solely by 
foreground emission. 

A two-component fit 

f v \ 20 

C‘ Q (e, V) = d° m (6) + C*<0) (-) ( 8 ) 

tests for the superposition of a CMB component with a single 
foreground component. Figure 4 shows the resulting decom- 
position for different combinations of WMAP frequency chan- 
nels. We detect foreground components in both IQ and IU, with 
best-fit spectral index fi = -2.4 ± 0.4 for IQ and fi = -3.0 ± 0.4 
for IU. It is likely that a single foreground component is re- 
sponsible for the signal in both IQ and IU; a joint fit yields 
fi = -2.7 ±0.3 with approximately equal amplitudes (at large 
angular scales) in curl- and curl-free modes, 0.5 ±0.1 jiK 2 at 
i/q = 41 GHz. The foreground detection does not depend solely 
on the 23 GHz data; we obtain similar (albeit noisier) results 


excluding the 23 GHz data from the fit. The fitted index is 
consistent with an origin from polarized synchrotron emission, 
and supports a picture in which we correlate polarized syn- 
chrotron emission with the total Galactic foreground, since the 
synchrotron, free-free, and dust morphologies trace each other 
to a great extent (Bennett et al. 2003c). It is thus not too sur- 
prising to find the best-fit foreground index in the temperature- 
polarization cross-correlation to be somewhat flatter than the 
spectral index for the synchrotron component alone. 

The fitted CMB component is robust against all combina- 
tions of frequency channels and fitting techniques. Note the 
agreement in Fig 4 between nearly independent data sets: the 
co-added QVW data (uncorrected for foreground emission) and 
the KKaQ data (corrected for foreground emission). We obtain 
additional confirmation by correlating /, Q , and U maps derived 
from a linear combination of the WMAP frequency channels 
designed to suppress foreground emission. We compute C IQ (9) 
using the “internal” linear combination 1 map of Bennett et al. 
(2003c), and similar minimum -variance linear combinations for 
the Q and U maps. The result is nearly identical to the CMB 
component derived from Eq. 8. 

3. POLARIZATION CROSS-POWER SPECTRA 

In a second analysis method, we compute the angular power 
spectrum of the temperature- polarization correlations using a 
quadratic estimator (Appendix A). The power spectrum is the 
Legendre transform of the two-point correlation function, and 
is more commonly encountered for theoretical predictions. We 
compute cJ E and cJ B individually for the each WMAP fre- 
quency band, using uniform weight for the temperature map 
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Fig. i . — Diagonal elements of the covariance matrix for the cJ E polarization cross-power spectrum. Points show the diagonal elements computed from 7500 
Monte Carlo simulations. The solid line shows the analytical model (Eq. 10). Note we multiply M// by [(/+ 1 )/( 27 t )] 2 to match the units in Figures 7 and 8. 


and noise weight for the polarization maps. We then combine 
the angular power spectra, using noise-weighted QVW data for 
/ > 21 where foregrounds are insignificant, and a fit to CMB 
plus foregrounds using all 5 frequency bands for / < 21. Since 
foreground contamination is weak, we gain additional sensitiv- 
ity in this analysis by using the Kp2 sky cut retaining 85% of 
the sky. 

We estimate the uncertainty in each l bin using the covariance 
matrix M for the polarization cross-power spectrum. Based 
on our analysis of the cf T covariance matrix (Hinshaw et al. 
2003b), the cJ E covariance matrix has the form along the diag- 
onal of 

Mu = < cf E cf E >-< cf E > 2 (9) 

_ (cf r + n T T/wi)(cf E +n EE /wi) + (cp) 2 
(2l + l)f sky f%, 

where htt and tiee are the TT and EE noise bias terms, w/ is the 
effective window function for the combined maps (Page et al. 
2003a), cJ T and cf E are the temperature and polarization angu- 
lar power spectra, = 0.85 is the fractional sky coverage for 
the Kp2 mask, and = f S ky/l ■ 14 for noise weighting. Figure 
5 compares the analytic expression for the diagonal elements 
of the covariance matrix to the mean derived from 7500 Monte 
Carlo simulations. The analytic form (Eq. 10) accurately de- 
scribes the simulations. We approximate the off-diagonal terms 
using the geometric mean of the covariance matrix terms for 
uniform and noise weighting (Hinshaw et al. 2003b), 10 

M;,/ ~ ( ) 0 ' 5 r A ; . (11) 

Figure 6 shows the off-diagonal terms r^i measured from 
Monte Carlo simulations. The largest contribution, -2.8%, is 
at A/ = 2 from the symmetry of our sky cut and noise coverage. 
The total anticorrelation is Sa/^o^/ = -0.124. Because of this 
anti-correlation, the error bars for the binned cJ E are slightly 
smaller than the naive estimate. A second method of estimating 
the errors relies on end-to-end simulations derived from sim- 
ulated rime-ordered data consisting solely of instrument noise 
(including the estimated contribution from 1// fluctuations). 


We have generated 1 1 “noise” sky maps each in Stokes /, Q , 
and U and compute the variance in TE directly from the vari- 
ance in the simulated signal. These two approaches yield errors 
that are consistent to better than 5%. 

Figure 7 shows the polarization cross-power spectra for the 
WMAP one-year data. The solid line shows the predicted signal 
for adiabatic CMB perturbations, based only on a fit to the mea- 
sured temperature angular power spectrum cf T (Spergel et al. 
2003; Hinshaw et al. 2003b). Two features are apparent. The 
TE data on degree angular scales (/ > 20) are in excellent agree- 
ment with a priori predictions of adiabatic models (Coulson 
et al. 1994). Other than the specification of adiabatic perturba- 
tions, there are no free parameters - the solid line is not a fit to 
cf E . The x 2 of 24.2 for 23 degrees of freedom indicates that 
the CMB anisotropy is dominated by adiabatic perturbations. 
On large angular scales (/ < 20) the data show excess power 
compared to adiabatic models, suggesting significant reioniza- 
tion. 

The WMAP detection of the acoustic structure in the TE 
spectrum confirms several basic elements of the standard 
paradigm. The amplitudes of the peak and anti-peak are a 
measure of the thickness of the decoupling surface, while the 
shape confirms the assumption that the primordial fluctua- 
tions are adiabatic. Adiabatic fluctuations predict a tempera- 
ture/polarization signal anticorrelated on large scales, with TE 
peaks and anti-peaks located midway between the temperature 
peaks Hu & Sugiyama (1994). The existence of TE correla- 
tions on degree angular scales also provides evidence for super- 
horizon temperature fluctuations at decoupling, as expected for 
inflationary models of cosmology (Peiris et al. 2003) 


4. TEMPLATE POWER SPECTRA 

Figure 7 demonstrates that the power spectrum of 
temperature-polarization correlations on degree angular scales 
can be predicted using the power spectrum of the temperature 
fluctuations alone. We use this for a third derivation of the 
TE cross-power spectrum, based on template matching in pixel 


10 Note that Hinshaw et al. (2003b) deft ne off-diagonal elements in terms of the inverse covariance matrix, which differs from /a/ by a sign. 




8 


Kogut et at. 


o 

o 


0.03 

0.02 p* 
0.01 
0.00 
-0.01 
-0.02 jr 


—0.03 


f- 

, ■ 1 1 ' , 

ll mI. 1 

' ] 


-40 


-20 


0 
A l 


20 


40 


FIG. 6. — Off-diagonal correlations r A / in the covariance matrix for the cJ E polarization cross-power spectrum, computed from simulations. All values are 
normalized to r&j = 1 at A/ = 0. The dotted line shows ta/ = 0 for comparison. The anti-correlation at A/ = 2 results from the spatial symmetry of the sky cut and 
noise coverage. 


space. On large angular scales, the signal-to-noise ratio for the 
temperature maps is much larger than one per multipole, while 
the S/N ratio for the polarization is much less than one. The 
likelihood function for the polarization measurement then has 
the simple form 

log C = (P-^a,^ red ) r AT 1 (P-^a ; ^ red ), (12) 

l l 

where P is the measured polarization signal (a 2 N pixe i vector), 
a/ = cJ E /cf T is the polarization fraction at each /, N is the pixel 
noise correlation matrix (a 2W pixe i x 2N pixe i matrix) and 

& pr Jn) = ^ 2 a , m ( 2 Y, m (n)+- 2 Y,m(n)) 

m 

Up Kd (n) = 2 Y lm (h)), (13) 

m 

Here ± 2 Yi m (h) are the spin harmonics, while a im are the mea- 
sured coefficients for an all-sky map of the CMB temperature. 
Imposing a cut to mask the Galactic plane introduces additional 
correlations; we avoid this by using the “internal” linear combi- 
nation temperature map (Bennett et al. 2003c) without impos- 
ing a sky cut. 

The maps Q pred and U pred represent the predicted polar- 
ization pattern based on the observed pattern of temperature 
anisotropy. We fit these template maps to the observed Q and 
U polarization maps to derive the polarization fraction a/ and 
thus the cJ E polarization cross-power spectrum. Minimizing 
the likelihood function yields the normal equations 



K w av =yi, 

(14) 

where 


yi = PN~r pre d 

(15) 

and 


Kll' = PpredN~ { Ppred- 

(16) 


These equations show the advantages of this approach. We 
compare the data with a template in pixel space, making it 
straightforward to include a spatially varying noise signal. We 
directly compare the measured polarization maps to a predic- 
tion based on the measured temperature maps, yielding a mea- 
surement of the TE polarization cross-power spectrum in the 


observed sky unaffected by cosmic variance. We can thus more 
easily compute the errors on the measured polarization fraction. 
The input temperature map (Stokes I) is already corrected for 
foreground emission (much simpler in pixel space where the 
unpolarized foregrounds are more easily measured), greatly re- 
ducing the foreground contribution to the cross-power spectra. 

We thus compute the temperature-polarization cross- 
correlation using three disparate techniques: the two-point an- 
gular correlation function, a quadratic estimator for the power 
spectrum in Fourier space, and a template fit in pixel space. 
All methods are in good agreement despite their very dissimi- 
lar treatment of the data. All methods show a significant excess 
of power for / <10. 

5. REIONIZATION 

WMAP detects statistically significant correlations between 
the CMB temperature and polarization. The signal on degree 
angular scales (/ > 20) agrees with the signal expected in adia- 
batic models based solely on the temperature power spectrum, 
without any additional free parameters. We also detect power 
on large angular scales (/ < 10) well in excess of the signal pre- 
dicted by the temperature power spectrum alone. This signal 
can not be explained by data processing, systematic errors, or 
foreground polarization, and has a frequency spectrum consis- 
tent with a cosmological origin. 

The signal on large angular scales has a natural interpreta- 
tion as the signature of early reionization. 11 Both the tempera- 
ture and temperature-polarization power spectra can be related 
to the power spectrum of the radiation field during scattering 
(Zaldarriaga 1997). Thomson scattering damps the temperature 
anisotropy and regenerates a polarized signal on scales compa- 
rable to the horizon. The existence of polarization on scales 
much larger than the acoustic horizon at decoupling implies sig- 
nificant scattering at more recent epochs. 

5.1. Reionization in a A CDM Universe 

If we assume that the ACDM model is the best description 
of the physics of the early universe, we can fit the observed 
temperature-polarization cross-power spectrum to derive the 


11 Although tensor modes can also generate TE correlations at large angular scales, tensor-to-scalar ratios r large enough to fi t the WMAP TE data are ruled out by 
the WMAP TT data (Spergel et al. 2003). 
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Fig. 7. — Polarization cross-power spectra c\ E for the WMAP one -year data. Note that we plot (/ + l)/27r cJ e and not /(/ + l)/27r cJ e . This choice emphasizes 
the oscillatory nature of cl E . For clarity, the dotted line shows q = 0. The solid line is the predicted signal based on the cf T power spectrum of temperature 
anisotropy - there are no free parameters. The TE correlation on degree angular scales (/ > 20) is in excellent agreement with the signal expected from adiabatic 
CMB perturbations. The excess power at low / indicates signifi cant reionization at large angular scales. 


optical depth r. We assume a step function for the ionization 
fraction x € and use the CMBFAST code (Seljak & Zaldarriaga 
1996) to predict the multipole moments as a function of optical 
depth. While this assumption is simplistic, our conclusions on 
optical depth are not very sensitive to details of the reionization 
history or the background cosmology. 

Figure 8 compares the polarization cross-power spectrum 
cJ E derived from the quadratic estimator to ACDM models 
with and without reionization. The rise in power for / < 10 
is clearly inconsistent with no reionization. We quantify this 
using a maximum-likelihood analysis 


£ oc 


exp(-4x 2 ) 

|M|‘/2 * 


(17) 


Figure 9 shows the relative likelihood £/Max(£) for the op- 
tical depth r assuming a ACDM cosmology, with all other 
parameters fixed at the values derived from the temperature 
power spectrum alone (Spergel et al. 2003). The likelihood 
for the 5-band data corrected for foreground emission peaks 
at r = 0.17 ± 0.03 (statistical error only): WMAP detects the 
signal from reionization at high statistical confidence. 

A full error analysis for r must account for systematic er- 
rors and foreground uncertainties. We propagate these effects 
by repeating the maximum likelihood analysis using different 
combinations of WMAP frequency bands and different system- 
atic error corrections. We correct C IQ (6) in each frequency 
band not for the best estimate of the systematic error templates, 
but rather the best estimate plus or minus one standard devia- 
tion. We then fit the mis-corrected C IQ (8 , z/) for a CMB piece 
plus a foreground piece (Eq. 8) and use the CMB piece in a 


maximum-likelihood analysis for r. The change in the best-fit 
value for r as we vary the systematic error corrections propa- 
gates the uncertainties in these corrections. Systematic errors 
have a negligible effect on the fitted optical depth; altering the 
systematic error corrections changes the best-fit values of r by 
less than 0.01. 

The largest non-random uncertainty is the foreground sepa- 
ration. We assess the uncertainty in the foreground separation 
by repeating the entire systematic error analysis (using both 
standard and altered systematic error corrections) with the fore- 
ground spectral index = -2.7 ± 0.3 shifted one standard de- 
viation up or down from the best-fit value. Table 2 shows the 
fitted optical depth r and goodness-of-fit statistic \ 2 f° r de- 
ferent data combinations and foreground spectral indices de- 
rived from the analysis of the two-point correlation function 
C IQ (0). The first set of rows shows values derived by simply 
co-adding the WMAP frequency channels, without any correc- 
tion for foregrounds. Data at 41, 61, and 94 GHz (Q, V, and W 
bands) where foregrounds are negligible show similar values 
for r; the x 2 ~ 63 for 57 degrees of freedom indicates that the 
data are in excellent agreement with reionized models. Adding 
additional low-frequency channels reduces the formal statisti- 
cal uncertainty, but introduces non-zero foreground contamina- 
tion as shown by the marked increase in x 2 . The next three 
sets of rows show the results when the data are separated into 
CMB and foreground components (Eq. 8). All data combina- 
tions are now in agreement; we obtain nearly identical values 
for r when fitting either the highest-frequency data set QVW or 
the lowest-frequency set KKaQ. The fitted optical depth is in- 
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Fig. 8. — WMAP Polarization cross-power spectra c T t E (fi lied circles) compared to ACDM models with and without reionization. The rise in power for / < 10 is 
consistent with reionization optical depth r = 0. 17 ±0.04. The error bars on WMAP data refect measurement errors only; adjacent points are slightly anti-correlated. 
The grey band shows the 68% confi dence interval from cosmic variance. The value at / = 7 is particularly sensitive to the foreground correction. 


sensitive to the spectral index: varying the spectral index from 
-2.4 to -3.0 changes the fitted values by 0.02 or less. We adopt 
r = 0.17 ±0.04 as the best estimate for the optical depth to 
reionization, where the error bar reflects a 68% confidence level 
interval including statistical, systematic, and foreground uncer- 
tainties. 

Spergel et al. (2003) include the TE data in a maximum- 
likelihood analysis combining WMAP data with other astro- 
nomical measurements. The resulting value, r = 0.17 ±0.06, 
is consistent with the value derived from the TE data alone. 
The larger uncertainty reflects the effect of simultaneously fit- 
ting multiple parameters. The TE analysis propagates fore- 
ground uncertainties by re-evaluating the likelihood using dif- 
ferent foreground spectral index. Since foreground affect only 
the lowest multipoles, the combined analysis propagates fore- 
ground uncertainty by doubling the statistical uncertainty in cJ E 
for 2 < / < 4 to account for this effect. 

5 .2. Model-Independent Estimate 

An alternative approach avoids assuming any cosmologi- 
cal model and uses the measured temperature angular cor- 
relation function to determine the radiation power spectrum 
at recombination. This approach assumes that the best es- 
timate of the three dimensional radiation power spectrum is 
the measured angular power spectrum rather than a model fit 
to the angular power spectrum. Given the observed tempera- 
ture power spectrum cJ T , we derive the predicted polarization 
cross-power spectrum cf E (§4), which we then fit to the ob- 
served TE spectrum as a function of optical depth r. We obtain 
r - 0A 6 ± 0.04, in excellent agreement with the value derived 


assuming a ACDM cosmology. We emphasize that the model- 
independent technique makes no assumptions about the cos- 
mology. The fact that it agrees well with the best-fit model from 
the combined temperature and polarization data (Spergel et al. 
2003) is an additional indication that the observed temperature- 
polarization correlations on large angular scales represent the 
imprint of physical conditions at reionization. The dependence 
on the underlying cosmology is small. 

5.3. Early Star Formation 

Reionization can also be expressed as a redshift z r assuming 
an ionization history. We consider two simple cases. For instan- 
taneous reionization with ionization fraction x e =\ at z< z r * the 
measured optical depth corresponds to redshift z r = 1 7 ± 3 . This 
conflicts with measurements of the Gunn-Peterson absorption 
trough in spectra of distant quasars, which show neutral hydro- 
gen present at z ~ 6 (Becker 2001; Djorgovski et al. 2001; Fan 
2002). Reionization clearly did not occur through a single rapid 
phase transition. However, since absorption spectra are sensi- 
tive to even small amounts of neutral hydrogen, models with 
partial ionization x e < 1 can have enough neutral column den- 
sity to produce the Gunn-Peterson trough while still providing 
free electrons to scatter CMB photons and produce large-scale 
polarization. Direct Gunn-Peterson observations only imply a 
neutral hydrogen fraction > 1% (Fan 2002). Accordingly, we 
modify the simplest model to add a second transition: a jump 
from x e = 0 to x e = 0.5 at redshift z r , followed by a second transi- 
tion from x e = 0.5 to x e = 1 at redshift z = 7. Fitting this model to 
the measured optical depth yields z r ps 20. In reality, reioniza- 
tion is more complicated than simple step transitions. Allow- 
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Fig. 9. — Likelihood function for optical depth r for a ACDM cosmology, using all 5 WMAP frequency bands fi tted to CMB plus foregrounds with foreground 
spectral index fi = -2.1. The best-fi t value is r = 0.17 with 95% confi dence range 0.09 < r < 0.28 after including systematic and foreground uncertainaties. 


ing for model uncertainty, the measured optical depth is consis- 
tent with reionization at redshift 1 1 < z r < 30, corresponding to 
times 100 < t r < 400 Myr after the Big Bang (95% confidence). 

Extrapolations of the observed ionizing flux to higher red- 
shift lead to predicted CMB optical depth between 0.04-0.08 
(Miralda-Escude 2002), lower than our best fit values. The 
measured optical depth thus implies additional sources of ion- 
izing flux at high redshift. An early generation of very massive 
(Pop DI) stars could provide the required additional heating. 
Tegmark (1997) estimate that 10~ 3 of all baryons should be in 
collapsed objects by z = 30. If these baryons form massive stars, 
they would reionize the universe. However, photons below the 
hydrogen ionization threshold will destroy molecular hydrogen 
(the principal vehicle for cooling in early stars), driving the ef- 
fective mass threshold for star formation to ^ 10 8 solar masses 
and im]oeding subsequent star formation (Haiman et al. 1997; 
Gnedin & Ostriker 1997; Tegmark 1997). X-ray heating and 
ionization (Venkatesan et al. 2001; Oh 2001) may provide a 
loophole to this argument by enhancing the formation of H 2 
molecules (Haiman et al. 2000). 

Cen (2003) provides a physically-motivated model of “dou- 
ble reic-nization” that resembles the two-step model above. A 
first generation of massive Pop III stars initially ionizes the in- 
tergalactic medium. The increased metallicity of the intergalac- 
tic medium then produces a transition to smaller Pop II stars, 
after which the reduced ionizing flux allows regeneration of a 
neutral hydrogen fraction. The ionization fraction remains at 
x e « 0.6 until the global star formation rate surpasses the re- 
combination rate at z - 6, restoring x e ~l. The predicted value 
r = 0.10 ± 0.03 should be increased somewhat to reflect the 
higher WMAP values for the baryon density Ct b and normal- 


ization erg (Spergel et al. 2003). 

6. CONCLUSIONS 

WMAP detects statistically significant correlations between 
the temperature and polarization maps. The correlations are 
inconsistent with instrument noise and are significantly larger 
than the upper limits established for potential systematic er- 
rors. The correlations are present in all WMAP frequency 
bands with similar amplitude from 23 to 94 GHz; fitting the 
data to a single power-law in frequency yields a spectral index 
fi - -0.4 ± 0.4, consistent with a CMB signal (ft = 0) and in- 
consistent with the measured spectral indices for Galactic fore- 
ground emission. A two-component fit to a superposition of 
CMB and Galactic foregrounds yields a positive foreground de- 
tection in both curl- and curl-free modes, with best-fit spectral 
index fi = -2.7 ± 0.3 consistent with synchrotron emission of 
amplitude 0.5 ± 0. 1 pK 2 antenna temperature at 4 1 GHz. 

The fitted CMB component is robust against different data 
combinations and fitting techniques. On small angular scales 
(0 < 5°), the WMAP data show the temperature-polarization 
expected from adiabatic perturbations in the temperature power 
spectrum. The data for £ > 20 agree well with the signal pre- 
dicted solely from the temperature power spectra, with no ad- 
ditional free parameters. 

The data show excess power on large angular scales (0 > 
10°) compared to the predictions based on the temperature 
power spectrum alone. The excess power is well described 
by early reionization at redshift z r = 20^°, corresponding to 
times t r = 1 801|q° Myr after the Big Bang (95% confidence). 
A model-independent fit to reionization optical depth yields re- 
sults consistent with the ACDM model. Our best estimate for 
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Table 2 

Reionization Optical Depth 2 


Data Set 

Method 

0 



f(>WMAP) Q 

VW 

Co-Add 


0-15^ 

63.0 

0.249 

QVW 

Co-Add 


0.18 ±0.03 

63.6 

0.233 

KaQVW 

Co-Add 


0.17 ±0.03 

91.7 

0.002 

KKaQVW 

Co-Add 


0.29 ±0.02 

352.1 

0.000 

KKaQ 

Co-Add 


0.28 ±0.01 

474.4 

0.000 

QVW 

Fit 

-2.4 

a 14+0.14 
u,h -0.04 

64.2 

0.218 

KaQVW 

Fit 

-2.4 

a 16+0.09 
u * 1D -0.03 

67.0 

0.162 

KKaQVW 

Fit 

-2.4 

0.18 ±0.03 

62.3 

0.265 

KKaQ 

Fit 

-2.4 

0.15 ±0.05 

60.6 

0.328 

QVW 

Fit 

-2.7 

a 14+0.13 
u * i 4 A ).04 

64.5 

0.214 

KaQVW 

Fit 

-2.7 

a 17+0.07 
u - i7 ^0.03 

66.5 

0.169 

KKaQVW 

Fit 

-2.7 

0.17 ±0.03 

68.1 

0.137 

KKaQ 

Fit 

-2.7 

0.14 ±0.05 

70.0 

0.112 

QVW 

Fit 

-3.0 

0 14+ 0 - 13 
ua ^-0.04 

64.5 

0.209 

KaQVW 

Fit 

-3.0 

0.18 ±0.05 

66.1 

0.172 

KKaQVW 

Fit 

-3.0 

0.17 ±0.03 

73.7 

0.061 

KKaQ 

Fit 

-3.0 

0.14 ±0.05 

79.1 

0.030 


"Optical depth r fi tted from d®($) for various combinations of data and foreground corrections in a ACDM cosmology. There are 57 degrees of freedom for each fi t. 
fc 68% confi dence statistical uncertainties 

^Fraction of 1000 simulations of reionized ACDM models with x 2 larger than WMAP value. 


the optical depth is r = 0.17 ±0.04 (68% confidence) where 
the error terms include statistical, systematic, and foreground 
uncertainties. This value is larger than expected given the de- 
tection of a Gunn-Peterson trough in the absorption spectra of 
distant quasars, and implies that the universe has a complex 
ionization history. 

The WMAP detection of early reionization opens a new fron- 
tier to explore the universe at redshift 6 < z < 30. WMAFs sen- 
sitivity to reionization is currently limited by instrument noise, 
both as direct statistical uncertainty and in the ability to better 
model and remove faint polarized foregrounds. Instrumental 
effects do not limit analysis of temperature-polarization cor- 


relations. The TE power spectrum and covariance matrix are 
available at http : / / lambda . gsf c . nasa . gov. We are 
currently performing a more complete set of systematic error 
analyses in the individual Q and U maps. A future data release 
will include full-sky polarization maps and polarization power 
spectra. 

The WMAP mission is made possible by the support of the 
Office of Space Sciences at NASA Headquarters and by the 
hard and capable work of scores of scientists, engineers, tech- 
nicians, machinists, data analysts, budget analysts, managers, 
administrative staff, and reviewers. 


APPENDIX 

QUADRATIC ESTIMATOR FOR TEMPERATURE-POLARIZATION POWER SPECTRUM 

We estimate the temperature-polarization power spectrum from pixelized sky maps using the following formalism. We begin by 
expanding the temperature and polarization fluctuations in generalized spherical harmonics 


T (») = a lmYlm(n) 

Itn 

(Al) 

Q(n) ± iU (n) = ^ a Tl ,im T 2 Yim(n) 

Iffl 

(A2) 

We then decompose the polarization fluctuations into E and B like pieces 


^±2 ,lm = Elm ± iBi m . 

(A3) 

We can use the basic properties of the spherical harmonics 


N Y lm = (-D (N+m) - N Y;_ m 

(A4) 

f MnhY^ih) = 5 Ui 

(A5) 
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to derive 

E lm = l -Jdn [ Q(n) ( ^Y'Kn) + - 2 Y,* m (h)) 

~>U(n){^Y;jn)-. 2 Yrjn)) ] 

B lm = ~\jdn [ U(n) ( ^Jn)+ ^(nj) 

~i Q(n ) ( T2 Y; m (n) - -iY,* m (n)) ] • (A6) 

We can now generalize the approach of Hivon et al. (2002) to estimate the coupling terms. We multiply the temperature and 
polarization maps by a weighting function 

7/m = Jdh w T (n) I(n)Y*,(n) (A7) 

Elm = { I dn w p (h) [Q(n)( 2 Y,' m (n)+ - 2 Y,* m (n)) 

-iU(n) ( 2 y,;(«)-_ 2 y, ,;(«))] (A8) 

Blm = dh w p (n) [U(n) ( 2 Y,* m (n)+ - 2 C(n)) 

-iQ(n)( 2 Y; m (n)~ - 2 Y? m (n))}. (A9) 

We expand the weighting function in spherical harmonics 


Mn) = y^w lm Yi m (n), 


Im 


and combine with equations Al - A3 to yield 


Tim = w[„ m „Ti< m i I dn Yii ml {h)Yin m n{h)Y* m (h) 

Vm'V'm" 


Eim= i E w /"»" [ dn Yi"„,"(n) ( 2 I'/'*'(n) 2 0«) + -2 Yi'm'(n) -2 C(n)) 

Z Vm'V'm" ^ J 

+iB,mjdh Yi" m "(n) ( 2 W(«) 2 ^* («)+ -2 W(«) - 2 ^(n)) 

B lm = - ^2 [^W f ^ 2C(") + -2^m'(«) -2^/m(«)) 

1 Vm'V'm" L ^ 

-iE lm j dn Y,n m n(n) ( 2 Y,, m ,(n) 2 Y? m (n)+ - 2 Y Vm dh) - 2 0«))J • 


We can then use 


I dfi N Y; m (h ) wY Vn ,(h) »"Yi„ m „(h) = (-\ ) N+m 


(2/+l)(2/'+l)(2/"+l) 


47 r 


1/2 


f / /' /" W / V /" \ 

\-NN'N") \-mrri m” ) 


to compute 


i TB 


¥ 

\ C? S ) 

After some algebra, the coupling terms reduce to 

.ttjt ( 2 /+ 1 )( 2 /'+ 1 ) 




C J 

c 




J 


r l L 

Jtb 


Jee 

CJ, 

\ c?, B } 


M, 


IV 


47T 




TEJE _ fofTBJB 


IV 


W 

(2/+ 1)(2/' + 1) 


87T 


E"H°oo) ! 


M EE,EE * *BBJIB 

w = M W 


(A 10) 


(All) 


(A 12) 


(A13) 


(A 14) 


(A 15) 
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(2/+l)(2/' + l)^ PP \f l /'/"W / Z'Z"\1 

[^-2 2 0 J + ^-2 2 0 J ] 


\6n 


i " 


\( i /'/"Y 

( 1 /' i"\] 

\-2 2 o y 

h \-2 2 0 J 


EEyBB _ M BB£E 

(2l+l)(2l 






Y / l'l"\ 

( l I'l" V 

(-2 2 0 J ' 

' (-2 2 0 J 


where 




(A 16) 

(A 17) 
(A 18) 


with a and b referring to either T or P. All of the other coupling terms are zero. Note that if we use different weighting functions for 
7\ Q and U, we increase the coupling between E and B modes. 

UNIFORM TEMPERATURE WEIGHTING 

If we use the full sky to compute the temperature spherical harmonic terms, then the cross-correlation term and its error matrix 
becomes particularly simple. For this case, = 1 / V 47T and all other coupling terms are 0. In this limit, the measured cJ E is just a 

constant times the true cJ E 


where 


The covariance matrix for these terms are diagonal. 


te _ c f E 
Cl ~ / 


cJ T cf E 




( 2 /+ 1)/ 2 

REFERENCES 


(Bl) 


(B2) 


(B3) 


Baccigalupi, C., Burigana, C., Perrotta, F., De Zotti, G., La Porta, L., Maino, 

D. , Maris, M., & Paladini, R. 2001, 4, 372, 8 
Barnes, C. et al. 2003, ApJ, submitted 
Becker, R. H. e. a. 2001, AJ, 122, 2850 

Bennett, C. L., Bay, M., Halpem, M., Hinshaw, G., Jackson, C., Jarosik, N., 
Kogut, A., Limon, M., Meyer, S. S., Page, L., Spergel, D. N., Tucker, G. S., 
Wilkinson, D. T., Wollack, E., & Wright, E. L. 2003a, ApJ, 583, 1 
Bennett, C. L., Halpem, M., Hinshaw, G., Jarosik, N., Kogut, A., Limon, M., 
Meyer, S. S., Page, L., Spergel, D. N., Tucker, G. S., Wollack, E., Wright, 

E. L., Barnes, C., Greason, M., Hill, R., Komatsu, E., Nolta, M., Odegard, 
N., Peiris, H., Verde, L., & Weiland, J. 2003b, ApJ, submitted 

Bennett, C. L. et al. 2003c, ApJ, submitted 
Bond, J. R. & Efstathiou, G. 1984, ApJ, 285, L45 

Brouw, W. N. & Spoelstra, T. A. T. 1976, Astron. Astrophys. Suppl. Ser., 26, 
129 

Bruscoli, M., Tucci, M., Natale, V., Carretti, E., Fabbri, R., Sbarra, C., & 
Cortiglioni, S. 2002, New Astronomy, 7, 171 
Cen, R. 2003, ApJ, submitted (astro-ph/02 10473) 

Coulson, D., Crittenden, R. G., & Turok, N. G. 1994, Phys. Rev. Lett., 73, 2390 

Djorgovski, S. G., Castro, S., Stem, D., & Mahabal, A. A. 2001, ApJ, 560, L5 

Fan, X. e. a. 2002, AJ, 123, 1247 

Gnedin, N. Y. & Ostriker, J. P. 1997, ApJ, 486, 581 

Haiman, Z., Abel, T., & Rees, M. J. 2000, ApJ, 534, 1 1 

Haiman, Z., Rees, M. J., & Loeb, A. 1997, ApJ, 484, 985 

Hinshaw, G. F. et al. 2003a, ApJ, submitted 

— . 2003b, ApJ, submitted 

Hivon, E., G6rski, K. M., Netterfield, C. B., Crill, B. P, Prunet, S., & Hansen, 

F. 2002, ApJ, 567, 2 


Hu, W. & Dodelson, S. 2002, ARA&A, 40, 171 
Hu, W. & Sugiyama, N. 1994, ApJ, 436, 456 
Hu, W. & White, M. 1997, Phys. Rev. D, 56, 596 
Jarosik, N. et al. 2003, ApJS, 145 
Kaiser, N. 1983, MNRAS, 202, 1169 

Kamionkowski, M., Kosowsky, A., & Stebbins, A. 1997, Phys. Rev. D, 55, 
7368 

Kosowsky, A. 1996, Annals of Physics, 246, 49 
Kovac, J. et al. 2002, ApJ, astro-ph/0209478 
Lubin, P. M. & Smoot, G. F. 1981, ApJ, 245, 1 
Miralda-Escude, J. 2002, ApJ, submitted (astro-ph/02 1 1071) 

Oh, S. P. 2001, ApJ, 553, 499 

Page, L. et al. 2003a, ApJ, submitted 

— . 2003b, ApJ, 585, in press 

Peiris, H. et al. 2003, ApJ, submitted 

Rees, M. J. 1968, ApJ, 153, LI 

Seljak, U. & Zaldarriaga, M. 1996, ApJ, 469, 437 

Spergel, D. N. et al. 2003, ApJ, submitted 

Tegmark, M., e. a. 1997, ApJ, 474, 1 

Tucci, M., Carretti, E., Cecchini, S., Nicastro, L., Fabbri, R., Gaensler, B. M., 
Dickey, J. M., & McClure-Griffi ths, N. M. 2002, ApJ, 579, 607 
Uyamker, B., Fuerst, E., Reich, W., Reich, P., & Wielebinski, R. 1998, Astron. 

Astrophys. Suppl. Ser., 132, 401 
— . 1999, Astron. Astrophys. Suppl. Ser., 138, 31 
Venkatesan, A., Giroux, M. L., & Shull, J. M. 2001, ApJ, 563, 1 
Zaldarriaga, M. 1997, Phys. Rev. D, 55, 1822 
Zaldarriaga, M. & Seljak, U. 1997, Phys. Rev. D, 55, 1830 



